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Abstract 



Standard statistical thermodynamic views of temperature fluctuations 
predict a magnitude (\/< (AT) 2 >/T) ps \/ (ks/C) for a system with heat 
capacity C. The extent to which low temperatures can be well defined is 
discussed for those systems which obey the thermodynamic third law in the 
form lim(T'_ > o) C = 0. Physical limits on the notion of very low temperatures 
are exhibited for simple systems. Application of these concepts to bound Bose 
condensed systems are explored, and the notion of bound Boson superfluidity 
is discussed in terms of the thermodynamic moment of inertia. 
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I. INTRODUCTION 



A problem of considerable importance in low temperature physics concerns physical limi- 
tations on how small a temperature can be well defined in the laboratory [1]. In what follows, 
we shall consider temperature fluctuations which define a system temperature "uncertainty" 

ST = .J< (T a - < T s >) > 2 , (1.1) 

whenever a finite system at temperature T s is in thermal contact with a large reservoir 
at bath temperature T. Only on the thermodynamic average do we expect the system 
temperature to be equal to the bath temperature; i.e. < T s >^ T. The fluctuation from 
this average result (quoted in the better textbooks discussing statistical thermodynamics^]), 
has the magnitude 

(£)- 

where C is the heat capacity of the finite system, and ks is Boltzman's constant. Since C 
is an extensive thermodynamic quantity, one expects the usual small fluctuation in temper- 
ature (ST/T) oc (l/x/TV) in the thermodynamic limit N — > oo, where N is the number of 
microscopic particles. However, in low temperature physics (for systems with finite values 
for N), temperature fluctuations are by no means required to be negligible. 
For those finite sized systems which obey the thermodynamic third law 

limC = 0, (1.3) 

one finds from Eqs.(1.2) and (1.3) that 

ST 

(— ) -> oo as T -> with iV < oo. (1.4) 

Eq.(1.4) sets the limits on what can be regarded as the ultimate lowest temperatures for 
finite thermodynamic systems; i.e. the temperature must at least obey ST «T. 

In Sec. II, the theoretical foundations for Eq.(1.4) will be discussed. In brief, the micro- 
canonical entropy of a thermodynamic system is given by 

S(E) = k B lnT(E), (1.5) 

where T(E) is the number of system quantum states with energy E. The micro-canonical 
entropy defines the system temperature T s via 




The thermal bath temperature T, which is not quite the system temperature T s , enters into 
the canonical free energy 

F(T) = -k B T lnZ(T), (1.7) 
where the partition function is defined as 

Z(T) = Tr(e- H ' kBT ) = £ T{E)e- E ' kBT . (1.8) 

E 
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The probability distribution for the energy of the system, when in contact with a thermal 
bath at temperature T, is given by 



N /Y(E) \ f-Ex ,F(T) - E + TS(E)s , . 

P(EiT > = W^W) = e *»( k B T - )• < L9) 

as dictated by Gibbs. Thus, the temperature T (of the thermal bath) does not fluctuate 
while system energy E does fluctuate according to the probability rule of Eq.(1.9). On the 
other hand, the system temperature T S (E) in Eq.(1.6) depends on the system energy and 
thereby fluctuates since E fluctuates. It is only for the energy E* of maximum probability 
M(ixeP(E;T) = P max — P{E*;T) that the system temperature is equal to the bath tem- 
perature T S (E*) = T. If the energy probability distribution is spread out at low thermal 
bath temperatures, then temperature fluctuations are well defined in the canonical ensemble 
of Gibbs. 

In Sec. Ill, temperature fluctuations are illustrated using the example of black body 
radiation in a cavity of volume V. For this case, it turns out that the thermal wave length 
A? of the radiation in the cavity, 

A^(^), (1.10) 

must be small on the scale of the cavity length V 1 ^ 3 ; i.e. A T « V 1 ^ 3 . For example, in a 
cavity of volume V ~ 1 cm 3 , the lowest temperature for the radiation within the cavity is 
of order T m j„ ~ 1 °K. It is of course possible to cool the conducting metal walls of a cavity 
with a length scale ~ 1 cm to well below 1 °K. However, this by no means implies that the 
radiation within the cavity can have a temperature well below 1 °K. The point is that at 
temperatures T s < 1 °K, there are perhaps only a few photons (in total) in the cavity. The 
total number of photons are far too few for the cavity radiation system temperature to be 
well defined. 

In Sec. IV, a confined system of Af atoms obeying ideal gas Bose statistics is discussed. 
Such systems can be Bose condensed, and are presently (perhaps) the lowest temperature 
systems available in laboratories. In the quasi-classical approximation, the free energy is 
computed in Sec.V. Questions concerning bounds on ultra-low temperatures are explored. 
Whether or not such Bose condensed atoms can exhibit superfluid behavior is discussed 
in Sec. VI. The superfluid and normal fluid contributions to the moment of inertia are com- 
puted. In the concluding Sec. VII, another simple system with fluctuation limits on ultra-low 
temperatures will be briefly discussed. 



II. THEORETICAL FOUNDATIONS 

Let 4>{E) denote some physical quantity which depends on the energy E of a physi- 
cal system. If the system is in contact with a thermal bath at temperature T, then the 
thermodynamic average may be calculated from 

<0>=5>(£;T)0(£), (2.1) 

E 
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where the probability P(E;T) has been defined in Eq.(1.9). Using the "summation by 
parts" [3,4] formula 

*E-^(P(E;T)<j>(Ej)=0, (2.2) 

E 



i.e. with a strongly peaked P(E; T) 



one finds 

-k B 



Employing Eqs.(1.6) and (1.9), 

Eqs.(2.4) and (2.5) imply the central result of this section 

-*b(%) = ((±-} f )A (2.6) 



dE' WT S T 



If we choose 4>(E) = 1, then Eq.(2.6) reads 



-) = -■ 

tj t 1 



(2.7) 



i.e. on average, the reciprocal of the system temperature is equal to the reciprocal of the 
thermal bath temperature. Thus, with fluctuations from the mean 

A0 = 0- < >, (2.8) 

» ( 1 \ 1 / 1 \ 1 1 / \ 

A (f ) = ¥.- %) = ¥.- r < 2 - 9 ' 

Eq.(2.6) reads 

If we choose in Eq.(2.10) the function to be 

-(i) = ^(§) = (5^) <*■»> 



(where C = (dE/dT s ) is the system heat capacity), then 

T^C 



A(I) 2 ) = (i i (^)). (2.12) 
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The standard Eqs.(l.l) and (1.2) follow from the more precise Eqs.(2.7) and (2.12) in the 
limit of small temperature fluctuations; i.e. 



7 rp2 

<(AT) 2 >^(^^) if ST = 7< (AT) 2 > « T. (2.13) 

The condition 5T « T is required in order that the canonical thermal bath tempera- 
ture be equivalent to the micro-canonical system temperature. If the micro-canonical and 
canonical temperatures are not equivalent, then the statistical thermodynamic definition of 
temperature would no longer be unambiguous. This raises fundamental questions as to the 
physical meaning of temperature. The view of this work is that in an ultra-low temperature 
limit, whereby 5T ~ T for sufficiently small T, the whole notion of system temperature is 
undefined, although the notion of a thermal bath temperature retains validity. 

III. BLACK BODY RADIATION EXAMPLE 

The heat capacity of black body radiation in a cavity of volume V with the walls of the 
cavity at temperature T is given by [5] 

C(Black Body) = fc B (^)(ii), (3.1) 

where A T is given by Eq.(l.lO). From Eqs.(1.2) and (3.1) it follows that the radiation 
temperature of a black body cavity of volume V is 



In order to achieve a well defined radiation temperature inside the cavity, 5T(Black Body) 
must be small on the scale of T or equivalently At << V 1 ^ 3 . As stated in Sec.l, this implies 
a minimum temperature of T min ~ 1 °K for a cavity of V ~ 1 cm 3 . 

IV. CONFINED IDEAL BOSE GAS 

The grand canonical free energy of an ideal Bose gas is determined by the trace [6] 

2 (T, n) = k B T tr ln(l - e 0^-h)/k B T^ ? ^ 
where h is the one Boson Hamiltonian and 

dE = -SdT - Nd^i (4.2) 
determines the number of Bosons Af. If the one Boson partition function is defined as 

q(T) = tr(e- h/kBT ), (4.3) 
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then the free energy obeys 

CO 1 rp 

E(T,n) = -k B Tj2(-)q( L )e nf * /kBT . (4.4) 

n=l 71 71 

The mean number of Bosons is 

OO rp 

X(T,/,)^^q(-)e n ^ T , (4.5) 

n=l 71 

and the statistical entropy is given by 

s(Tfi) = _ ^(r,^) + ^(r,^) ) + fer g ( J_ )? , ( r )e „^ Ti (46) 

where g'(T) = {dg(T)/dT}. 

Of considerable theoretical [7,8] experimental[9,10,ll] interest is the bound Boson in an 
anisotropic oscillator potential, 



where 




Eqs.(4.3) and (4.7) imply 

"( r )-n( 1 _ e -U. r )- («) 

The heat capacity may be defined by 

which must be calculated numerically. 

Shown in Fig.l is a plot of the heat capacity (in units of Nks) versus temperature (in 
units of the critical temperature T c ) for the case of M = 2 x 10 3 and M = 2 x 10 6 particles. We 
choose, for experimental interest[12], the frequency eigenvalues (ui/2n) = (uo 2 /2t{) = 320Hz, 
and (cu 3 /2n) = ISHz. 
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FIGURES 



Heat Capacity vs. Temperature 



N=2,000 (Dotted) N=2,000,000 (Solid) 
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FIG. 1. The heat capacity (in units of NUb) is plotted as a function of temperature (in units 
of T c ) for M = 2 x 10 3 atoms (dotted curve) and M = 2 x 10 6 atoms (solid curve). 

Figure 1. 

For finite A/", there is strictly speaking no Bose-Einstein condensation phase transition. 
The critical temperature T c is therefore defined as that temperature for which the heat 
capacity reaches the maximum value C max = C(T C ). Although phase transitions are defined 
in mathematics only in the thermodynamic limit Af — > oo, for all practical purposes, a 
quasi-classical approximation of Eq.(4.1) in the form 

S(T, n) = k B T J J (j^|)Zn(l - e (M-h(p,r))/k B T^ ^ quasi _ dassica i^ ( 4 . 10 ) 

where 

h(p,r) = ^+ 1 -Mr.u 2 .r, (4.11) 

does yield a Bose-Einstein condensation phase transition whose heat capacity is sufficiently 
accurate for values of M ~ 10 5 or higher. Thus, we regard recent experiments on Bose atoms 
confined in a magnetic bottle to be probing a physical Bose-Einstein ordered phase. Let us 
consider Eq.(4.10) in more detail. 
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V. QUASI-CLASSICAL BOSE CONDENSATION 

In order to evaluate Eq.(4.10) we employ the quasi-classical form [13, 14] of Eqs.(4.3) and 
(4.11); i.e. 

where u> = (uoi^uo?,) 1 ^ . From Eqs.(4.4) and (5.1), it follows that Eq.(4.10) evaluates to 



71=1 

From Eqs.(4.2) and (5.2), the number of particles obeys 



"< T 'ri = (^)' ( " <0) ' (5 - 3) 

With the usual definition of the ^-function 

OO 1 

CW = Eb). ne(s)>l, (5.4) 



n=l - n * 



the Bose-Einstein condensation critical temperature is 



The non-zero number of Bosons (below the critical temperature) in the condensate state is 
given by 

A/- (T)=A/"{l-(^) 3 }, (T<T C ). (5.6) 
Finally, the entropy below the critical temperature 

S( r)=« flC (4)(^) 3 = «B(^)(|) 3 , (T<T C ), (5.7) 

obeys the thermodynamic third law limp^o) $(T) = 0. The heat capacity in the Bose- 
Einstein Condensed phase is then given by 

Employing Eqs.(1.2) and (5.8) we find that the temperature uncertainty below the critical 
temperature obeys 

Thus, for A/" ~ 10 5 one may safely consider the temperature of the ordered phase to be well 
defined in the range T c > T > T min , where T min ~ 0.05 T c . The open question as to whether 
the ordered phase is a superfluid may now be considered. 
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VI. SUPERFLUID FRACTION OF THE BOUND BOSON SYSTEM 



The notion of a superfluid fraction in an experimental Bose fluid (such as liquid He 4 ) 
may be viewed in the following manner: Suppose that we pour the liquid into a very slowly 
rotating vessel and close it off from the environment. The walls of the vessel are at a bath 
temperature T, and the vessel itself rigidly rotates at a very small angular velocity Q. In 
the "two-fluid" model[15,16], the normal part of the fluid rotates with a rigid body angular 
velocity Q, which is the same as the angular velocity of the vessel. On the other hand, the 
superfluid part of the fluid does not rotate. The superfluid exhibits virtually zero angular 
momentum for sufficiently small f2. The total fluid moment of inertia tensor / is defined by 
the fluid angular momentum L = I-fl (as S7 — > 0). We here take the limit Q — > 0, to avoid 
questions concerning the effects of vortex singularities on the superfluid. The normal fluid, 
which rotates along with the rotating vessel, contributes to the fluid moment of inertia. The 
superfluid, which does not rotate with the vessel, does not contribute to the moment of 
inertia. Thus, the geometric moment of inertia, 



where p(r) is the mean mass density of the fluid (at rest), overestimates the physical moment 
of inertia eigenvalues when the fluid is actually a superfluid. The normal fluid contributes 
to the moment of inertia and the superfluid does not do so in the limit fi — * 0. Below, we 
consider in detail the moment of inertia of the bound Bose gas. 

For the bound Bose system, we consider a mesoscopic rotational state [17] with a thermal 
angular velocity Q. The rotational version of Eq.(4.2) reads 




(6.1) 



dE n = -SdT - Afdfi - L-dil, 



(6.2) 



where L is the bound Boson angular momentum. Eq.(4.11) gets replaced by 



Mp,r) 




+ -Mr-LJ 2 -r - n-(r x p), 



(6.3) 



so that Eq.(5.1) now reads 




■ b t _ 



(q(T)/y/Det(l 



(6.4) 



where the matrix Co 2 is written in Eq.(4.8) and 



n 2 



( (nl + n§) 



(ft? + n 2 ) -n 2 n 3 

-fl 2 Q 3 {n 2 + Q 2 ) ) 



(6.5) 



From Eqs.(4.4) and (6.4), it follows that 




(6.6) 
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The fluid moment of inertia tensor has the matrix elements 



dL d 2 ^ 



(6.7) 



Eqs.(4.8), (6.5), (6.6) and (6.7) imply (in the unordered phase) 



, (T>T C ). 



(6.8) 



In the unordered phase, obeying Eq.(5.2), one finds that Eq.(6.8) is precisely what would 
be expected from a normal fluid with geometric moment of inertia 



4 = J d 3 rp(r)(r% - r^), (T > T c ), 



(6.9) 



where p(r) is the mean mass density of the atoms. 

In the ordered phase (T < T c ), the moment of inertia of the particles over and above the 
condensate is given by Eq.(6.8) with jj, = 0, i.e. 



rexcitation 



\UJ 2 U). 



V 



2 ^3 






\ 

(i + i)/ 



(T < T c ). (6.10) 



The question of superfluidity concerns the magnitude of the moment of inertia of those 
particles within the condensate. For T < T c , we use the notation that I denotes the mo- 
ment of inertia of the excited Bosons, and J represents the moment of inertia of the Bose 
condensate. If the moment of inertia of the particles in the condensate were zero, then the 
condensate particles would all be "superfluid" . 

Let ipo(r) be the normalized (/ d 3 r\ip (r)\ 2 = 1) Bose condensation state. From the 
geometric viewpoint, the moment of inertia of the condensate would be given by 



geometric = ^ J fo^f^ _ ^ 



i.e. 



jgeometric 



hM 



/(- + -) o 

o + 










(6.11) 



(6.12) 



(- + -) 

V^l ^2/ 

The physical Bose condensate moment of inertia tensor is in reality 

jphysical =N Q J2 ( - t^M^ >< j^Mh. > + < j^Mjl >< t^Mh. > \ (6 13) 

where 1 = —ih(r x V). One may derive Eq.(6.13) by treating the rotational coupling Ah = 
— f2 • 1 to second order perturbation theory in the energy Ae (f2) as Q — > 0. Eq.(6.13) 
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evaluates to 



jphysical _ 



hAf 



\102<JJ'3 (0)2+0)3) 



V 



















(6.14) 



V 0)10)2(^1+^2) / / 



If there exists an axis of symmetry, then the condensate moment of inertia corresponding 
to that axis vanishes. For an experimental example, if lo\ = 002 7^ ^3, then the 3— axis is an 
axis of symmetry and jP^ ysical = as implied by Eq.(6.14). In this sense, for the ideal Bose 
gas, the superfTuid fraction is the same as the condensate fraction in Eq.(5.6), 



T 3 

1] super fluid{T) = |l - ( — ) }. 



(6.15) 



If we were to employ an axis which is not a symmetry axis, i.e. if the angular momentum 
about that axis were not conserved, then even the superfluid would give some contribution 
to the moment of inertia. 

Finally, including the scattering length will have some effect on the magnitude of the 
superfluid fraction[18]. 



VII. CONCLUSIONS 

We have employed the heat capacity C in the relationship 



(t) * m) <"> 

in order to place lower bounds on possible system temperatures. The minimum system 
temperatures were estimated based on the notion that the temperature is "well defined" 
only if temperature fluctuations are small 5T «T. The third law of thermodynamics, in 
the form lim( T ^ ) C = 0, dictates that the condition 5T « T is harder to achieve as the 
temperature is lowered. 

For example, one may achieve low temperatures by adiabatic demagnetization[19]. The 
thermodynamics of the method is well illustrated by a system of M non-interacting two level 
particles, each with possible energies 

E± = ±A. (7.2) 

Such a physical system might consist of M nuclear one half spins in a magnetic field. The 
mean energy of such a system, 

E = -ATA tanh( A/k B T) , (7.3) 
implies a heat capacity C = (dE/dT) given by 

C-Afk ( {A/kBT) Y (7 4) 
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The temperature fluctuation of a system of non-interacting two level particles follows from 
Eqs.(7.1) and (7.4) to be 

( 5T \~( 1 u cosh(A/fc B T) ^ 

M W (v^)( (A/k B T) )• (7 - 5) 

Thus, for a system of A/" ~ 10 6 two level particles, it is only possible to achieve a 
temperature as low as T min ~ 0.2 (A/k B ). For a system of A/" ~ 10 6 trapped Bosons, as 
shown in this work, it is possible to achieve a temperature as low as T min ~ 0.05 T c . It is 
evident on the grounds of temperature fluctuations alone that the bound boson systems are 
the more likely example for the very lowest temperatures. Indeed, this has turned out to be 
true in the laboratory. 
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